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Abstract 

I-V characteristics of the high T^ superconductor Bi2Sr2CaiC208 shows a strong 
hysteresis, producing many branches. The origin of hysteresis jumps is studied by 
use of the model of multi-layered Josephson junctions proposed by one of the authors 
(T. K.). The charging effect at superconducting layers produces a coupling between 
the next nearest neighbor phase-differences, which determines the structure of hys- 
teresis branches. It will be shown that a solution of phase motions is understood as 
a combination of rotating and oscillating phase-differences, and that, at points of 
hysteresis jumps, there occurs a change in the number of rotating phase-differences. 
Effects of dissipation are analyzed. The dissipation in insulating layers works to 
damp the phase motion itself, while the dissipation in superconducting layers works 
to damp relative motions of phase-differences. Their effects to hysteresis jumps are 
discussed. 
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1 Introduction 

I-V characteristics of the high Tc superconductor Bi2Sr2CaiC208 (BSCCO) shows a strong 
hysteresis, producing many branches [jl], 0, ^. Many experiments have been reported to 
investigate properties of hysteresis branches [^, ^ |, |^, ||, |], |I^, and it has turned out 
that there is rich physics in properties of I-V characteristics. Substructures are found in 
hysteresis branches, and are discussed in relation with phonon effects |^. The d-wave 
effect is also discussed in higher voltage region, where the quasi-particle tunneling current 



becomes important P, |T0 



The present paper is concerned to the origin of hysteresis jumps forming the main 
branches. Instead of attributing the origin of hysteresis branches to inhomogeneity of 
critical currents, there are several theoretical attempts to regard its origin as an intrinsic 
property of strong anisotropy in high T^ superconductors. High T^ superconductors have 
a layered crystal structure. Specially, BSCCO has a strong anisotropy. Then a model of 
multi-layered Josephson junctions, which has an alternating stack of superconducting and 
insulating layers, is considered as a proper model to describe electronic properties of those 
systems. One of main differences from the case of a single junction is the appearance of an 
interlayer coupling among phase-differences. Since the thickness of the superconducting 
layer is of the order of 3-6A, influence among adjacent junctions is not negligible. 

To explain the main branch structure, the following mechanisms are proposed. A 
mechanism of the interlayer coupling induced by the charging effect at the superconducting 



layer has been proposed by one of the author (T. K.) [|TI], and a branch structure of I-V 
characteristics is obtained |T^. It is pointed out in ref. |jTB|| that a nonequilibrium effect 
in superconducting layers also mediates an inter-layer coupling, and I-V characteristics 
are discussed |1^, |[. When the magnetic field is applied, the inductive coupling becomes 
important, and its effects have been discussed [^, |15|, [16], |1^, p!8[| . 

In this paper, we investigate in detail the scheme of hysteresis jumps appearing in 
main branches without magnetic filed. Since the inductive coupling among adjacent junc- 
tions does not exist, the charging effect at superconducting layers becomes important. A 
charging at a superconducting layer induces the change of electric fields in the neighbor- 
ing insulating layers and induces an inter-layer coupling, since the time-derivative of the 
phase-difference is related with the voltage. We will show that the origin of the multi- 
branch structure is the nonlinear effect among phase- differences through the charging 
effect. Various disturbances, such as dissipation in insulating layers, dissipation in su- 
perconducting layers, noise, a surface proximity effect affect schemes of hysteresis jumps. 
We will investigate how the schemes of hysteresis jumps are affected by the above men- 
tioned disturbances. The dissipation in the superconducting layer will be introduced 
phenomenologically through the charge relaxation. Schemes of hysteresis jumps are ob- 
tained by increasing applied current gradually up to certain value and then decreasing. 
Obtained equations show that the resistive dissipation in the insulating layer works to 
damp the motion of phase-difference in each junction, while the charge relaxation in the 
superconducting layer works to damp relative motions of phase differences among neigh- 
boring junctions. This fact affects the way of formation of hysteresis jumps in the process 
of increasing and decreasing applied current. It will turn out that the dissipation in su- 
perconducting layer and noise enhance the occurrence of jumps in both increasing and 
decreasing current, although the main branch structure is determined by the coupling 



through the charging effect. By comparing the result of numerical simulation with the 
experiment of a small stack [0, we will show validity of the present model. 

In the next section, the necessary formula for calculation is summarized. The dissi- 
pation in superconducting layers is introduced phenomenologically. In Sect. 3, results 
of numerical simulation are presented. The origin of hysteresis jumps is investigated. It 
will be shown that, at a position of a jump, a stationary solution changes its form. We 
present a comparison with the experiment of ref. and show validity of the charging 
mechanism. Sect. 4 is devoted the conclusion. 



2 Josephson Junction Stacks 

In this section we summarize the multi- Josephson junction model of high T^ supercon- 



ductors |Tl|]. Let us consider the A^ -|- 1 superconducting layers, numbered from to A^. 
We denote the gauge invariant phase difference of the {I — l)-th and l-th superconducting 
layer by ip{l), and voltage by V{1). The widths of the insulating and superconducting 
layers are denoted by D and s, respectively. At the edges, the effective width of the 
superconducting layer may be extended due to the proximity effect into attached lead 
metals. The widths of the 0-th and A^-th superconducting layers are denoted by sq and 
Sjv, respectively. 

We assume that physical quantities are spatially homogeneous on each layer. This 
assumption is applicable in cases of no applied magnetic field and small sample size in the 



a-b direction. According to the procedure presented in ref. |TT|, we have the following 
equation for the total current J flowing through each junction as 

Jc Vo ^pOt Vo 

The current J is normalized by the critical current Jc, ic{l) = Jd})/ Jc with Jc(/) being 
the critical current for the /-th junction. The time t is normalized by the inverse of the 



plasma frequency, ujp = J^^^^^^, with e being dielectric constant of the insulating layer. 

The voltage V{1) is normalized by Vq = -^. The parameter /3 is given by /3 = j^. In 
this paper we assume 

jcil) = 1 for all / , (2.2) 

for simplicity. The relation between time-derivative of phases and voltages at junctions 
is given by 

(2.3) 
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A= -a l + 2a -a ■■■ , (2.4) 

where 

(2.5) 



with /i^ being given by the relation between the charge density and scalar potential Aq as 

m 

47rp = -/i-2(Ao + ^^0), (2.6) 

where is the phase of the order parameter defined on the superconducting layer. In 
the present analysis, we assume that the current fiows only along the c-axis, say, the 
z-direction. Since the current jz in superconducting layers is given by 
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the current conservation law leads to 
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The gauge condition ( p.8| ) is called phason gauge |]T9[, and is determined in such a way 
that the response of the current and charge to the space-time variation of phase appears 
in a gauge invariant form. It is shown in ref. |1^ that a nonequilibrium effect caused 
by the tunneling current induces the dissipation effect in /i~^. Since the superconducting 
layer is thin, the current also receives effects from boundaries between superconducting 
and insulating layers, which may also works as dissipation to the superconducting current. 
When the effect of the d-wave superconductivity is considered, one expects also certain 
dissipation due to the presence of the gap-less region. Such effects have not been much 
studied yet. Then we introduce phenomenologically the dissipation in superconducting 
layer in a Lorentz form as 
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This corresponds to the following phenomenological form of charge relaxation. 
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By considering the expression a in Eq. ( |2.5| ), we replace a by 

a -^ a{l — iujr) 
with a being given by Eq. ( p.5|) replaced iJ? by /ig- From Eqs. 
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where A is the matrix A in Eq. 
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As the above derivation shows, the charging effect at superconducting layers is the 
main mechanism inducing the inter-layer coupling through the electromagnetic interac- 
tion, when the spatial homogeneity in the a-b plane is satisfied. A charge at a supercon- 
ducting layer modifies electric fields in the neighboring insulating layers, and the change 
of electric fields affects the time-derivative of phase differences. Phonon effect considered 
in ref. ^ may be introduced through the frequency dependence of the dielectric constant 
e, which is not considered in this paper. The equation ( p.l3| ) means that the dissipation 
in the insulating layer works as a damping of phase motions at each junction, and that the 
dissipation in the superconducting layer works as a damping of the relative phase motion. 
In the next section, we present results obtained by numerical simulation of Eq. ( 2.13|) . 



3 Results of Numerical Simulation 



We have performed numerical simulations of Eq. ( 2.13| ) by use of the 5-th order of the 



predictor corrector method. In all calculation, we choose the parameter a = 1.0, P = 0.2. 
The time step dt is chosen as Updt = 1.0 x 10~^. 

The parameter j3 is related to the conductivity of the insulating layer. Although 
it works to lead a system into a stationary state and determines the slope of the I-V 
characteristics, its effect to a scheme of hysteresis jumps is not sensitive in this parameter 
range. The parameter a plays an important role to determine hysteresis branches due to 
non-linear effects among phase-differences. The value of a in BSCCO was estimated in 



ref. |TT[] as 1 < a < 3 and, by use of e = 25 for La2_xSra,Cu04, the value of a = 2.7 was 
obtained. There is a discussion that e for BSCCO may be smaller, and from a microscopic 
theory, a is estimated much smaller value as 0.2 [Q. Also it was reported that the equal 
spacing of I-V branches observed in experiments of BSCCO is reproducible, with a smaller 
choice of a < 1 [0. Here we use a typical value a = 1.0. 

By use of Eqs. ( |2.3| ) and (^.l]), the relation for obtaining a voltage for each junction 
is written as 

ujp or ^, Vq Jc 

The average of voltage V{1) is given by 

_ I f J- max 

V{1) = / dtV{l) . (3.2) 
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The total DC-voltage V is obtained by 

N 

V = Y.V{1). (3.3) 

1=1 

In Fig. 1, I-V characteristics is presented for A^ = 10, sq/ s = sn/s = 2, rup = 0.1. 
The current J is gradually increased with the step dJ/Jc = 0.01 up to J/Jc = 2.0, and 
then gradually decreased. The I-V curve shows the jump at J/Jc = 1-0 and, after two 
jumps, it linearly increases up to J/Jc = 2.0. In the current decreasing process, there 
appear three jumps when the current becomes smaller than J/Jc = 0.4. 



In Fig. 2, we show the branch structure in the I-V characteristics corresponding to 
Fig. 1. This is obtained by decreasing the current when a jump occurs in the current 
increasing process, and by increasing the current when a jump occurs in the current 
decreasing process. It should be noted that the positions of jumps labeled by 1 and 5 
lie on the same branches, while the positions 2 and 4 are close but lie on the different 
branches. 

In order to understand this branch structure, we show, in Fig. 3, the voltage distribu- 
tions on each junction just after the jumps. Fig. 3(a) is for the current increasing process 
and Fig. 3(b) is for the current decreasing process. For the line labeled 1, two junctions at 
the edges have higher voltage. The solution shows that the phases with higher voltages in- 
creases approximately linearly in time, corresponding to the phase rotating motion, while 
other junctions have oscillating phase motions. In fact, when the time average of -^ is 
plotted in Fig. 4, junctions with rotating and oscillating phases are clearly identified. As 
can be seen in Fig. 3, the number of junctions with rotating phase increases as 2, 4, and 
6 due to the symmetric arrangement. The number of phase rotating junctions is two for 
the voltage lines labeled by 1 and 5, four for 2 and 4. The lines 1 and 5 shows the same 
pattern of voltage distribution and therefore form the same branch in Fig.3, while the line 
2 and 4 show the different pattern, though the number of the phase rotating junction is 
same. 

With this analysis, we see that the hysteresis jumps are associated with the change 
of the distribution pattern of rotating phase motions. If the l-th junction has a rotating 
phase, the time average of -^ is constant and that of sin(y9(/)) is zero. If the l-th. 
junction has an oscillating phase, the time average of -^ is zero and that of sin((y9(/)) is 
constant. Let us consider that Nr{> 0) junctions have rotating phases. Considering the 
above facts, we can obtain the equation to determine the volt age- current relation from 
the time-average of Eqs. (|2.1|) and ( p.l3| ) as 



P^ = N^-J2 Uk) < sin ^(/o) > , (3.4) 
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where < ■ • ■ > indicates time average, /q and Iq are the labels of junctions with oscillating 

phases, 

s s 

a(l) = 1 + a — , a{N) = 1 + a — , otherwise a{l) = 1 . (3.6) 

So Sn 

As is shown in Fig. 3, even if the number of phase-rotating junctions is same, slightly 
different branches are formed according to the pattern of distribution of rotating phases. 
That is, branches are grouped by the number of phase-rotating junctions, Nr, and different 
patterns of distribution of phase-rotating junctions lead shghtly shifted branches. The 
Eqs. ( p.4|) and ( p.5|) are solved in a form 

P^ = n^. (3.7) 

In Fig. 5, we present n versus Nji, the number of rotating phases. The small dots are for 
all possible patterns. The edge junctions mostly rotate, and rotating phases have mostly 



oscillating phases in the next neighbors. Taking into account these facts, we have the 
larger circles with the restrictions, 1) pattern is symmetric, 2) two edges have rotating 
phases and 3) there are no patterns with more than three successive rotating phases. 
From the figure we can see that branches appears roughly with equal spacing but that 
those at higher voltage distribute denser. In the numerical simulation, it seems that only 
special patterns are realized. The simulated result in Fig. 2 gives n = 3.21 for N^i = 2, 
n = 6.21 and 6.47 for Nji = 4, and n = 8.59 for N^ = 6. What stability conditions should 
be imposed in addition to Eqs. ( |3.4} ) and ( ^.5[ ) is still an open question. 

In order to see how the scheme of hysteresis jump changes with various conditions, 
we show, in the following, structures of branches in I-V characteristics. In these analyses, 
the current is increased up to J/Jc = 2 and decreased to zero. Branches are picked up 
when a voltage jump occurs. Fig. 6 is for finite junction stack A^ = 10, and Sq = sj\f = 2.0 
is chosen. The parameters are a = 1.0, /? = 0.2. In Fig. 6, (a) is for UpT = 0.0 without 
noise, (b) is for UpT = 0.1 without noise, (c) is for UpT = 0.0 with noise, and (d) is for 
UpT = 0.1 with noise. The noise is introduced through the small randomness of the time 
derivative of the phase at each time step, S{^-^) = 1-0 x 10~^. Fig. 7 is for N=10 
with the periodic boundary condition. Fig. 7(a) is for UpT = 0.0 without noise, (b) is for 
UpT = 0.1 without noise, (c) is for UpT = 0.0 with noise, and (d) is for UpT = 0.1 with 
noise. 

When a noise is introduced, possible patterns are not necessarily symmetric, so that 
the number of branches increased. In addition, an instability to cause a hysteresis jump 
is much easily enhanced. The effect of r also enhance hysteresis jumps, which may be 
understood from the fact that it works as damping of the relative motion. Specially, it 
enhances hysteresis jumps in the current increasing process. 

In the periodic cases, this tendency is more enhanced as can be seen from comparison 
with Fig. 7(a) and others. In a long stack, where surface effects are less, we can say 
that the dissipation effect in the superconducting layer or noise play an important role to 
induce hysteresis jumps. 

In Fig. 8, we present the comparison with our numerical result and the experiment 
of ref. [^. We choose N = 8 from their arrangement of the sample. In order to explain 
large extension of hysteresis branches in the current increasing process, we assume that 
the critical currents Jc are reduced at the edges. We choose Jc = OATmeV, since the 
numerical analysis shows that the Jc is situated roughly at the middle of the second 
branch. From the data, we have Jc{0)/Jc = 0.4, Jc{N)/ J^ = 0.9 and other Jc{l)/ Jc = 1-0. 
The parameter a is chosen 1.0, by assuming s = 6A, D = 6A, fx =2Aand e = 10. In ref. 
2^, it is reported that the longitudinal plasma energy is of order of O.bmeV. Then Vq is 



0.25'meV. The P is estimated from the branch 4; that is, the branch 4 has four rotating 
phases and is on the line passing J = 0.9mA, V = 60meV. In Eq. (|3.7| ) we choose n = Q 
from the result of Fig. 5, and P is estimated as /? = 0.05. Since we do not introduce any 
other parameters, we set sq/s = 1.0, sn/s = 1.0. The parameter UpT is chosen as 0.5, 
which is obtained by trial and error to produce as many branches as possible. The inset 
is from ref. [^. The labels of branches correspond to those in the inset. The numbers in 
the brackets indicate the number of rotating phases. The values of J, at which hysteresis 
jumps occur, fit very well with experimental values of the inset. Hysteresis jumps are 
induced when certain instability develops in the non-linear equation of ( p.l3|) . Though 



each junction, except for one at edges, has the same Jc, a critical J for each branch is 



different and is controlled mainly by a. The agreement of this critical J- values observed 
in Fig. 8 supports the validity of the present charging mechanism. In Fig. 8 certain 
branches in the increasing current are missing. How hysteresis jumps occurs depends of 
effects of disturbances. There may be certain noises to make development of instability 
easier. The length of branches depends on choice of parameter. In decreasing current, 
there are a few discrepancies. In numerical calculation, the I-V curve is straight, while 
the experiment shows a curvature in the low current region. In addition, the positions 
of jumps are different. Such discrepancy may relate effects of d-wave superconductivity, 
such as the presence of the quasi-particle tunneling, for example, but it is beyond the 
discussion of the present work. 

4 Conclusion 

In this paper we have studied the I-V characteristic branches in high Tc superconductors 
by the multi-layered Josephson junction model. We have shown that the charging effect 
in superconducting layers causes an inter-layer coupling among phase-differences and that 
hysteresis jumps are induced as a result of non-linear effect among phase-differences. The 
appearance of hysteresis branches is intrinsic in high Tc superconductors. On this aspect, 
we note the experimental report that, even in a sample with negligible Jc- inhomogeneity, 
a branch structure is observed [|10| . 

Due to the weak interlayer coupling arising from the charging effect, equations for 
phase-differences form a non-linear coupled differential equation. Its solution is classified 
by the number of rotating phases. The change of its number occurs intrinsically in this 
nonlinear equation, which leads to formation of I-V characteristic branches. Although the 
number of rotating phases is same, different patterns of distribution of rotating phases 
lead slightly shifted branches. It is still an open question and a future problem which 
patterns are more easily realized among many possible patterns. 

In this system, the resistive dissipation in the insulating layer works as damping of 
phase motion in each junction and the dissipation caused by charge relaxation in super- 
conducting layers works as damping of relative phase motions. The effects of dissipation, 
noise, and surface proximity to structure of hysteresis jumps in I-V characteristics are 
investigated. The dissipation in superconducting layer and noise enhance the occurrence 
of jumps in both increasing and decreasing current. 

We have presented the comparison with the experiment in Fig. 8. The agreement of 
the critical J for each hysteresis branch shows validity of the charging mechanism as the 
origin of the inter-layer coupling. 

Finally, we comment on effects of the d-wave superconductivity. The fact that the 
Josephson current is given by the phase-difference is not modified in the d-wave super- 
conductivity with a single component. Therefore, the main conclusions in the present 
paper are not modified. As was pointed out in the text, the d-wave superconductivity 
gives, due to the presence of a gapless region, an additional contribution to dissipation 
in superconducting layers and additional quasi-particle tunneling current even in the low 
voltage region. Effects of such contributions may be interesting future problems. 
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Figure captions 

Fig.l I-V characteristics 

Parameters are a = 1.0, /3 = 0.2, N = 10, f = ^ = 2.0 and LUpT = 0.1. 

Fig. 2 I-V characteristics and branch structure 

Parameters are a = 1.0, /3 = 0.2, AT = lo, ^ = ^ = 2.0 and cUpT = 0.1. 

Fig. 3 Voltage distribution 

Parameters are a = 1.0, /3 = 0.2, iV = 10, f = ^ = 2.0 and ujpT = 0.1. (a) is 
for increasing current, and (b) is for decreasing current. 

Fig. 4 Average -^-^ distribution 

Parameters are a = 1.0, (3 = 0.2, iV = 10, ^ = ^ = 2.0 and UpT = 0.1. (a) is 
for increasing current, and (b) is for decreasing current. 

Fig. 5 Slope n vs. the number of rotating phases N^ 

Slope n is defined by (3y- = Uj-- Parameters are a = 1.0, A^ = 10, ^ = ^ = 
2.0. Small dots correspond to all possible patterns, while larger circles are with 
restrictions presented in the text. 

Fig. 6 I-V characteristics and branch structure for a finite stack 

Parameters are a = 1.0, /3 = 0.2, N = 10, f = ^ = 2.0. {&)ujpT = 0.0 without 

noise, (h)u!pT = 0.1 without noise, {c)ujpT = 0.0 with noise, and {d)ujpT = 0.1 

with noise. 
Fig. 7 I-V characteristics and branch structure for a periodic stack 

Parameters are a = 1.0, (3 = 0.2, A^ = 10, ^ = ^ = 2.0. (a)w„r = 0.0 without 



s 



noise, {h)LUpT = 0.1 without noise, (c)cc;pr = 0.0 with noise, and {d)ujpT = 0.1 
with noise. 
Fig. 8 I-V characteristics 

Inset is from ref . . The labels of the branches are identified with those of the 
inset. The numbers in the bracket are the numbers of rotating-phases. 
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